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METHOD FOR CHARACTERIZING AN OPTICAL FIBER LINK 

Related Applications 

[0001] This is a continuation application under 35 U.S.C. § 120 of WO 
2005/041449 Al, filed as PCT/BE2004/000154 on October 27, 2004, which is hereby 
incorporated by reference. This application claims the benefit of the filing date of U.S. 
Provisional Patent Application No. 60/514,625, filed October 27, 2003, which is hereby 
incorporated by reference. 

Field of the invention 

[0002] The present invention is related to a method for characterizing an 
optical fiber link by analysis of a Rayleigh backscatter signal, from which the spatial 
distribution of the beat length, the coupling length and the PMD (Polarization Mode 
Dispersion) along the link can be determined. 

State of the art 

[0003] The use of dispersion-shifted and dispersion compensating fibers has 
minimised the effect of chromatic dispersion on the bandwidth of an optical link. 
Polarization mode dispersion (PMD) has therefore become the most serious limiting 
factor in high-speed optical communication systems. There exist several measurement 
techniques of PMD, but they only allow for measuring the global value and do not give 
information about the distribution of PMD along the fiber length. Measuring the spatial 
distribution of PMD is, however, important. Such a measurement would indeed allow to 
locate the bad trunks in an optical link, causing high PMD value. This can be essential in 
the frame of a network maintenance and upgrade. 

[0004] The PMD of a fiber depends on two parameters: the beat length (Lb) 
and the coupling length (L c ). The beat length depends on the birefringence of the fiber 
and the coupling length is related to the phenomenon of polarization mode coupling. 
Measuring the spatial distribution of PMD therefore involves measuring the distribution 
of both the beat and coupling lengths. 

[0005| Most of the techniques proposed for measurement of the distributed 
PMD are based on polarization-optical time domain reflectometry (POTDR). The concept 
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of POTDR was introduced by A. Rogers in the 80's (see A. Rogers, 'Polarization-optical 
time domain rejlectometry: a technique for the measurement of field distributions', 
Applied Optics, vol. 20, pp. 1060-1074, 1981). It basically involves measuring the 
polarization properties of the Rayleigh backscatter signal when an optical pulse 
propagates down the fiber. Some POTDR-based techniques have since been developed 
for PMD characterization. Most of them allow a measurement of the beat length 
distribution along an optical link and, therefore, do not allow the complete determination 
of the PMD. However, Sunnerud presented a POTDR-based set-up enabling the 
measurement of the accumulation of PMD along the fiber (s ee H. Sunnerud et al, 
'Polarization-Mode Dispersion measurements along installed optical fibers using gated 
backscattered light and a polarimeter', J. of Lightwave Technol, vol. 18, pp. 897-904, 

2000) . This method requires POTDR measurements at several wavelengths. A technique 
based on the measurement of the degree of polarization by using a POTDR set-up has 
also been presented in 'Distributed detection of high-PMD sections on installed fibers 
using a polarization-OTDR\ M. Leblanc, Proc. OFMC'01, pp. 155-162, 2001, and is 
related to the patent application US 2003/0174312. This technique does not allow 
quantification of the PMD, but enables to obtain distributed information about the level of 
PMD. Moreover this technique requires the complete measurement of the state of 
polarization of the backscattered signal. More recently, Galtarossa described another 
technique for the PMD measurement for which the coupling length is determined from 
the correlation of the birefringence vector (see 'Measurement of birefringence correlation 
length in long, single-mode fibers', Galtarossa et al., Opt. Lett., vol. 26, pp. 962-964, 

2001) , which also requires the measurement of the polarization states of the backscattered 
signal. Moreover, this technique requires a quite small duration of the pulses launched in 
the optical fiber: 5-10 ns. 

[0006] In 'Polarization mode dispersion mapping in optical fibers with a 
Polarization-OTDR', M. Wuilpart, G. Ravet, P. MegretandM. Blondel, Photon. Technol. 
Lett., vol 14, pp. 1716-1718, 2002, the authors describe an analysis of the POTDR signal 
which allows a mapping of the PMD along an optical fiber link by quantifying the PMD 
on each fiber of the link. The main advantage of that technique is that it does not require 
the complete measurement of the backscattered polarization state evolution and only uses 
a linear polarizer at the fiber input, which is quite simple to implement. The determination 
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of the beat length, the coupling length and finally the PMD is based on the analysis of the 
statistical properties of the extrema present in a POTDR trace. 

[0007] In said paper by Wuilpart et al. the spatial distribution of the linear 
birefringence 5 is characterized by a Rayleigh statistical distribution along the fiber length 
and it is assumed that the rate of change of the birefringence angle q is described by a 
Gaussian distribution with zero mean and standard deviation a depending on the 
polarization mode coupling strength. The circular birefringence is neglected. The 
corresponding measurement set-up is shown in Fig. 1. The OTDR 110 pulses (10 ns) 
externally modulate a 1553 nm distributed feedback (DFB) laser 120 via a pulse 
generator 130 and an electro-optic modulator (EOM) 140. After amplification, the pulses 
are launched into the fiber through an acousto-optic modulator (AOM) 150, which 
suppresses the amplified spontaneous emission noise of the erbium doped fiber amplifier 
between two successive pulses. A linear polarizer 160 is placed at the fiber input and a 
polarization controller (CP) 170 is used to obtain the maximum power after the polarizer. 

[0008] Fig. 2 shows three simulated POTDR traces corresponding to the 
POTDR arrangement and obtained for different values of the beat length and an identical 
a value. One can easily observe that the number of maxima strongly depends on Lb. 

[0009] The calculation of the beat length is based on the measurement of the 
number of maxima of the POTDR trace. Fig. 3, which has been obtained by simulation, 
presents the number n of maxima per unit of length in function of the beat length for 
different values of a. This figure shows that n varies not only with L B but also with the 
strength of the polarization mode coupling (a). In order to estimate Lb from n, the curve 
corresponding to a = 10 (the dashed curve in Fig. 3) is used. This choice indeed 
minimises the maximum possible error if one assumes that a is smaller than 15 
degrees/m. Consequently, this method allows only to give an approximate value of the 
beat length Lb. A relationship between the coupling length L c and a can be derived and 
one can obtain: 

0.879 

L c~ — — 
a 

[0010] A statistical approach is then used for determining the coupling length 
from the POTDR trace. Fig. 4 shows different POTDR traces that were obtained by 
simulation for a beat length of 25 m and for different values of a: 2, 4 and 8 degrees/m, 
respectively. One can observe that the lower envelope of the curves varies in a smoother 
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way when a decreases and therefore when Lc increases. Defining APj as the absolute 
value of the difference in ordinate between two successive minima of the POTDR trace P 
one can write : 

where Zj and z ]+ \ correspond to the distances of the i th and the i+l lh minimum, 
respectively. 

[0011] A new parameter £ can then be calculated as 

<f =< A/> >, 

where <.>j denotes the mean value for i varying from 1 to the number of minima of the 
POTDR trace. % is a measure of the speed variation of the lower envelope of the 
backscattered signal along the fiber length. Fig. 5 (obtained by simulation) shows the 
evolution of £ according to a for different values of the beat length Lb. The values are 
here expressed in dB because % has been calculated directly from the POTDR trace. 

[0012] The procedure for the determination of Lc is the following. From the 
POTDR trace, the number of maxima and each APj are measured. An approximate value 
of the beat length is then deduced by the method previously described and £, is calculated 
by averaging APj. Finally a can be deduced from § by using the curve of Fig. 5 
corresponding to the correct beat length and then Lc is calculated by means of the relation 
Lc=0.879/a 2 stated above. Because L B is an approximate value, the calculated coupling 
length Lc will also be an approximation. 

[0013] From the measurement of the beat length and the coupling length, it is 
easy to determine an approximate value of the PMD value by 
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where c is the light velocity in a vacuum, X the measurement wavelength and L the fiber 
length. It should be emphasised that the final PMD value is approximate: this is due to the 
approximation done on the beat length measurement. This process can be repeated for 
each contribution in the POTDR signal corresponding to the different fibers of the optical 
link: it is therefore possible to map the PMD along this link. 

[0014] The method as described above has some important drawbacks and 
limitations. It only yields approximate values of L B > L c and the PMD. Further, the 
presented POTDR traces were assumed ideal. It supposes ideal characteristics of the 
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various components of the experimental set-up. In practice, the set-up is characterized by 
some imperfections, which cause a distortion of the ideal POTDR signal. Therefore the 
method as such cannot be directly applied. The POTDR signal is also affected by noise: 
detector noise and residual coherence noise, which adds a series of minima and maxima 
to the POTDR trace. The final value of the measured PMD can therefore be erroneous. 

Summary of Certain Embodiments 

[0015] Certain embodiments provide a method for analysing a backscatter 
signal for determining in an accurate way the spatial distribution of the beat length, the 
coupling length and the polarization mode dispersion. 

[0016] Certain embodiments relate to a method of characterizing an optical 
fiber link by its beat length, coupling length and polarization mode dispersion distribution 
comprising 

- sending a pulsed signal along the optical fiber link and measuring the backscattered 
signal (preferably a POTDR signal), after passing through a polarizer, 

- deriving the length of the optical fiber, the average power difference between two 
successive minima of the backscattered signal and the number of maxima per unit 
length, 

in an iterative way determining a beat length interval and an interval for the 
polarization mode coupling parameter, until the length of the intervals is below a 
predetermined value, yielding a value for the beat length and the coupling length, and 

- calculating the polarization mode dispersion. 

[0017] Preferably the POTDR signal is an ideal POTDR signal. 

[0018] Advantageously the POTDR signal is the convolution of an ideal 
POTDR signal and a signal depending on the pulse shape. In an alternative embodiment, 
the POTDR signal further is convoluted with a signal taking into account the time jitter 
introduced by the measurement set-up. 

[0019] In a specific embodiment a smoothing algorithm is applied to the 
POTDR signal. 

[0020] Another embodiment relates to a method for characterizing an optical 
link comprising a concatenation of several fibers, wherein the method as previously 
described is applied to each fiber. 
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[0021] Another embodiment relates to the use of a method as described above 
to locate the position of polarization mode dispersion sources within an optical fiber link. 

[0022] In an advantageous embodiment, the method as described is used in 
telecommunication networks. Alternatively, it can be applied in fiber sensing 
applications. 

Short description of the drawings 

[0023] Fig. 1 is a block diagram of an example set-up for measuring a POTDR 

trace. 

[0024] Fig. 2 is a diagram of three simulated POTDR traces for a given a and 

for different L B values. 

[0025] Fig. 3 is a diagram of the number of maxima per unit length versus L B 

for different values of a. 

[0026] Fig. 4 is a diagram of simulated traces for a given beat length and 
different values of a. 

[0027] Fig. 5 is a diagram of the evolution of £ as a function of a for different 
values of Lb. 

[0028] Fig. 6 is a diagram of the determination of the beat length interval. 

[0029] Fig. 7 is a diagram of the determination of the interval for a. 

[0030] Fig. 8 is a diagram of the distortion of the POTDR signal due to the 
imperfections of the set-up (pulse shape + time jitter). 

[0031] Fig. 9 is a diagram of the evolution of n with the beat length when the 
imperfections of the set-up are taken into account. 

[0032] Fig. 10 is a diagram of the evolution of £ with a when the 
imperfections of the set-up are taken into account. 

[0033] Fig. 1 1 is a diagram of the simulated and measured probability density 
function of A. 

[0034] Fig. 12 is a diagram representing the window algorithm that can be 
used for calculated the local PMD. 

Detailed description of certain embodiments 

[0035] It is first shown how the polarization mode dispersion from a POTDR 
signal can be determined in an accurate way (Lb is not approximately determined 
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anymore) by analysing the statistical properties of its extrema. Suppose first that a falls 
within the range [0,15] degrees/m, which means that the coupling length is assumed to be 
larger than 12.8 m. This a interval has been chosen such that it corresponds to the values 
measured on real fibers. 

[0036] From a POTDR signal, the fiber length L, £ and the number of maxima 
per unit of length, n m , are easily measurable. After measuring n m , it is possible to 



determine a first interval 



B. 



low up 



for the beat length of the fiber under test so that: 



1} 1} 

R 5 R 

"low "up 



L l R and L are determined by the intersections of the curves n(L B ,a) corresponding 

"low up " ' 

to a = 0 and a = 1 5 degrees/m with the straight line n = n m as shown on Fig. 6. 

[0037] When the first beat length interval is determined, the curves £>(cj,Lb) 

can be used to define an interval o } 0 w > a up for the a value of the fiber under test so that: 

a e [a hwi a up 



low 
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<?iow and & U p are determined by the intersections of the curves ^(a,L B ) for L B ~~ • g 

and Lb = L with the straight line £ = where £ m is the £ value of the fiber under 
test measured from the POTDR signal. This principle is illustrated in Fig. 7. 

[0038] After the determination of crj ow and &l p , the curves n(LB,a) for 

G = a bw and & = (J u P can be calculated. A new beat length interval is then deduced as 
previously explained by Fig. 6, but using the curves obtained for <r = <j} ow and o ' = <J X up 



instead of a = 0 and a = 15 degrees/m. 



[0039] This new beat length interval denoted 



b 



B. 



law tJ up 



is included in 



[} i} 

"low "up 



and gives a smaller range for the beat length of the fiber under test. Fig. 7 is 



2 2 

then used again to calculate a new a interval p/mv»°"wpj by using the curves £(a,L B ) 



obtained for 
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L B - L D and l b - 



up 



This new interval gives a smaller range for the a value of the fiber under test. This 
process is repeated until small enough intervals are obtained for the beat length Lb and for 
a. 

[0040] When the beat length L B and a are determined, the PMD is finally 
calculated by means of the relationship: 



[0041] If the optical link to be characterized comprises a concatenation of 
several fibers, the method is applied to each fiber of the link and their beat and coupling 
lengths are therefore determined independently of each other. In this way, the PMD 
mapping is possible along the link. 



several advantages : 

it requires only the access to one end of the optical fiber link to be characterized; 
it is based on a POTDR arrangement, but does not require the complete measurement 
of the complete state of polarization of the backscattered signal. Only a simple 
polarizer is required; 

it involves a POTDR measurement at only one wavelength; 

it gives accurate (not approximate) values for the beat length, the coupling length and 



[0043] In practice, the final measured signal is also affected by noise. This 
noise is due to the detector noise and to the presence of a residual coherence noise. As the 
determination of the PMD is based on the detection of extrema, it is obvious that this 
noise will affect the results. 

[0044] As already mentioned before, the measured signal is also affected by 
the imperfections of the experimental set-up like the finite pulse shape and the time jitter 
of the electronic equipment. 

[0045] These phenomena (noise and imperfections of the measurement set-up) 
lead to a distortion of the POTDR signal. Therefore it cannot be used directly to deduce 
the PMD map. 




[0042] 



The method for analysing the backscatter signal as described offers 



the PMD. 
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[0046] The width of the optical pulses launched into the fiber under test is 
finite and can be non-rectangular. Consequently, the backscattered power measured at the 
POTDR corresponds to an average of the ideal signal over a small spatial window 
depending on the pulse shape. The resulting signal measured at the POTDR is in fact the 
convolution product between the ideal signal and the shape of the pulse: 

where P pu ise(z) is the power measured by the POTDR according to a certain pulse shape. 
P(z) is the ideal POTDR signal and P S ha P e(z) is the power distribution along the pulse. 
Hence the Fourier transform of P pu i se (z) becomes: 

P pulse 

where P pu ise(k), P(k) and P S hape(k) are the Fourier transforms of P pu ise(z), P(z) and P S ha P e(z), 
respectively. 

[0047] Moreover the electronic equipment is not perfect: it is subject to a time 
jitter, e.g., the time delay between two successive pulses is not constant. This 
phenomenon strongly affects the averaging process of the POTDR. If we suppose that the 
number of pulses involved in the averaging process of the OTDR is large enough and that 
the time jitter is a random variable characterized by a Gaussian distribution of zero mean 
and standard deviation a t j, the signal detected at the POTDR can be written as: 

*■ 

P m {z) = P pulse {z)®e 2a « 
where P m (z) is the power measured by the POTDR taking into account the pulse shape 

and the time jitter. Hence the Fourier transform of P m (z) becomes: 

P m (k) = P pulse {k)o tj e 2 
The pulse shape and the time jitter have for effect to distort the measured backscattered 
signal. An example is shown on Fig. 8, which has been obtained for parameters of the 
experimental set-up. The signals (a) and (b) correspond to the ideal and affected signals, 
respectively. One can clearly observe that the signal is strongly distorted. 

[0048] Because of these imperfections, the method described above cannot be 
used directly: the curves n(L B ,a) and £>(a,L B ) are indeed based on ideal POTDR signals 
obtained by simulations. Therefore the effects of the pulse shape and the time jitter should 
be included in the calculation of n(L B ,a) and £(a,L B ). In the case of the parameters of the 
experimental set-up and a pulse duration of 10 ns, these new curves are shown on Fig. 9 
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and 10. In order to deduce these curves, ideal POTDR traces have been simulated and the 
effects of the imperfections have been taken into account by using the above expressions 
for Ppuhe(z) and P m (z). Consequently, when the imperfections of the experimental set-up 
are taken into account, the PMD map of an optical link can still be determined. The same 
reasoning as the one explained previously has to be applied on the POTDR signal by 
using the new set of curves n(L B> a) and £(a,L B ) described in Figs. 9 and 10. As the pulse 
shape is now taken into account, it is possible to launch a larger pulse in the fiber (50 ns, 
for example), which will increase the system dynamics. The maximum measurable length 
will therefore be greater. 

[0049] As already indicated, the final measured signal is also affected by 
noise. This noise is due to the detector noise and the presence of a residual coherence 
noise. Its main effect is to add a series of minima and maxima to the POTDR trace. The 
values of n and £ cannot thus be correctly measured and the final calculated PMD will be 
erroneous. In order to solve this problem, a smoothing algorithm is applied to the 
backscattered POTDR signal P m (z). This algorithm comprises generating the signal P s (z) 
such that: 



P s (z) is therefore the mean value of P m (z) on a distance interval L s around z. The number 
of maxima and the distances Zk corresponding to the locations of the minima of the 
POTDR trace are then calculated from the smoothed signal in order to ignore the 
influence of the small peaks generated by the noise. Afterwards, when the positions Zk are 
known, £ is calculated from P m (z). When noise is present, the number of extrema of P s (z) 
depends of the smoothing parameter L s . So, the main problem is to choose the correct 
smoothing parameter to be applied to the POTDR signal. 

[0050] The location of the minima of the smoothed signal can be used to 
estimate the probability density function of A, denoting the difference in z between two 
successive minima. This estimation can then be compared to the density probability 
function of A of a POTDR trace obtained by means of simulations taking into account the 
parameters L B and a calculated from the smoothed signal. If the probability density 
functions of A obtained from P s (z) and from the simulation do not match, it means that 
the noise contribution is significant. The smoothing parameter L s will therefore be chosen 
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such that the corresponding probability density function obtained from measurements can 
be correctly fitted by the probability density function derived from simulations: L s is 
therefore chosen such that the noise contribution is reduced. Fig. 1 1 illustrates this 
concept. Figures 11(a), 11(b) and 11(c) show the probability density functions obtained 
by measurement (circles) and by simulation (crosses) for a given optical fiber and for an 
L s value of 3, 4 and 5 m, respectively. One can clearly observe in this case L s = 3 m gives 
the best fitting. 

10051] Certain embodiments provide a technique which allows to locate more 
accurately the position of the PMD source within the fiber. The basic idea involves 
applying the described method on a distance window (1 km, for example), which slides 
along the fiber length. By comparing the obtained PMD's on several successive windows 
as shown on Fig. 12, it is possible to locate more accurately high PMD sources within the 
optical fiber. 

[0052] POTDR can also be used in fiber sensing. The polarization properties 
indeed depend on several external effects like strain, temperature, electric and magnetic 
fields. Strain, twist, temperature, electric and magnetic fields affect the birefringence 
within the fiber and consequently, the mean beat length. Certain embodiments, which 
measure the spatial distribution of the beat length, therefore allow to measure the spatial 
distribution of these external effects after calibrating the sensing fiber. 
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